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Relativistic Brownian motion can be inexpensively demonstrated on a graphene chip. The inter- 
play of stochastic and relativistic dynamics, governing the transport of charge carrier in graphene, 
induces noise-controlled effects such as (i) a stochastic effective mass, detectable as a suppression of 
the particle mobility with increasing the temperature; (ii) a transverse ratchet effect, measurable as 
a net current orthogonal to an ac drive on an asymmetric substrate, and (iii) a chaotic stochastic 
resonance. Such properties can be of practical applications in the emerging graphene technology. 
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As astroparticle data often require a relativistic anal- 
ysis, the challenge of a consistent formulation of rela- 
tivistic statistical thermodynamics is so much intriguing 
as timely [ij. Such a challenge encompasses the notion 
of relativistic Brownian motion as well. In particu- 
lar, issues like the relativistic generalization of inherently 
non-local thermodynamic quantities, such as heat and 
work, or the very concept of an equilibrium heat bath, 
have not been unambiguously settled, yet [H The 
state of the art in this field thus calls for experiments ca- 
pable to assess and validate the existing theoretical ap- 
proaches. This presently constitutes a nearly impossible 
task, due to the scarcity of the available relativistic data 
sets, which often also lack the necessary accurateness. 
Typically, such data either are being obtained from cos- 
mic rays Q or could be, in the future, from expensive 
high-energy experiments. To overcome this difficulty, 
we propose here an alternative and affordable route to 
demonstrate the physics of relativistic Brownian motion 
under controllable laboratory conditions, namely on a 
suitable graphene chip. 

Graphene is the only artificially crafted material 
known so far, which is truly two-dimensional (2D) Q. 
It is essentially a monolayer of carbon atoms packed in a 
honeycomb lattice and isolated from the bulk. The fact 
that graphene has a 2D structure makes its electronic 
properties rather unique. Charge carriers propagating 
through such a lattice are known to behave as relativistic 
massless Dirac fermions 0, @] . On the quantum mechani- 
cal level, they are described by the 2D analog of the Dirac 
equation with the Fermi velocity, vp = 10^ cm/s, replac- 
ing the speed of light in the Dirac description. Close to 
the so-called Dirac points, the energy-momentum rela- 
tion of the elementary excitations, often associated with 
quasi-particles, is linear: the velocity of the charge car- 
riers is thus always collinear with the momentum, its 
modulus being constant and equal to vp. 

Our proposal relies on the observation that for strong 
applied electric fields or at high enough tempera- 
tures, quantum transport effects, otherwise dominant in 



graphene, are effectively suppressed and the motion of 
charge carriers can be satisfactorily described by classical 
(i.e., non-quantum) relativistic equations. Such a semi- 
classical description has been invoked 0, for instance, 
to successfully interpret graphene's unusually broad cy- 
clotron resonance 0] • 

As detailed with this Letter, the combination of 
stochastic and relativistic dynamical effects governing the 
carrier motion in graphene results in some unexpected 
transport phenomena. The source of stochasticity is pro- 
vided by the finite temperature, whereby the charge car- 
riers are being scattered by phonons, lattice defects and 
buckling, and the sample boundaries, which leads to an 
equilibrium redistribution of their energies. We empha- 
size, however, that a noisy environment in the chip can 
also be generated by external noise sources, controllable 
or not, independent of the temperature, e.g., in the form 
of current or voltage fluctuations. Moreover, in graphene 
the Fermi velocity, which plays the role of the light speed 
in the relativistic equations for its charge carriers, is rela- 
tively low; one can then analyze driven relativistic Brow- 
nian motion in the presence of electromagnetic fields by 
ignoring otherwise hardly tractable relativistic retarda- 
tion effects. Under these conditions, driven relativistic 
Brownian particles become accessible to both theory and 
experiment [jl. [lot. 

Model ~ The scheme of a simple graphene chip is 
sketched in Fig.[lja). An undoped graphene sheet sits 
on a periodic sequence of parallel electrodes with alter- 
nate constant potentials, respectively, and U~ , with 
> U~. By tuning the distance between electrode 
pairs within a fixed spatial period, L, the ensuing stag- 
gered electric potential, U{x), can be modulated at will, 
symmetrically or asymmetrically, with amplitude U in 
the X direction, while being identically zero in the y 
direction. An example of the directed potential U{x) 
is drawn in Fig.[lja). Additionally, the charge carri- 
ers are driven by a spatially homogeneous electric field, 
E{t) = {Ex{t),Ey{t)), with sinusoidal components Ex{t) 
and Ey{t) acting along the x and the y axis, respectively. 
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FIG. 1: (Color online) Panel (a): Sketch of a graphene chip 
(carbon atoms arranged in a honeycomb lattice) with a set of 
several electrodes placed parallel to the y-axes. The elec- 
trodes, kept at a constant potential, U"^ (green) and U~ 
(blue), create an asymmetric electric potential, U{x), shown 
as a solid line. External ac electric fields, Ey{t) and Ex{t), 
drive the charge carriers in the x and y directions, respec- 
tively. Panel (b): Strength of the spectral peak, S, for a 
relativistic particle moving on a uniform substrate, U{x) = 0, 
parallel to the ac drive E = (0, A cos {uyt)): S/Sq vs. T/To, 
with So = (A/'yujy)^ and To = v\'^l2ujy (see text), for dif- 
ferent LOy (legend). Simulation parameters are: A = 20, and 
= 7 = fc = 1. The dashed line represents the asymptotic 
decay law, T^^, discussed in the text. 



The random dynamics of the charge carriers in a 
graphene sheet can be described by a set of coupled 
Langevin equations for the components of the 2D mo- 
mentum, p = {px,Py), subject to the relativistic dis- 
persion relation e — vpPo, where po = (p^ + Py)^ , £ 
is the particle energy and vf the Fermi velocity. The 
form of such phenomenological Langevin equations is 
determined by the condition that the chosen particle- 
reservoir coupling must lead to the same equilibrium 
momentum distribution as predicted by the fully mi- 
croscopic (Hamiltonian) theory. At low temperatures, 
the equilibrium distribution function of the quasiparti- 
cle (electron) energy obeys the Fermi-Dirac statistics, 
p{£) = l/(cxp[e/fcT] -1-1), where T is the temperature 
and the chemical potential was set to zero to indicate 



that the number of the quasiparticles is undetermined. 
However, the operation of the graphene chips considered 
here typically requires carrier energies e ^ C/, so that 
for sufficiently high energies, kT ^ e ^ U, p{£) is con- 
veniently approximated by the relativistic Jiittner distri- 
bution, p(j>) - exp(-£(^/fcr) 0, [iH. A viable set of 
Langevin equations proven to be consistent with a 2D 
Jiittner distribution reads |2| 



Px = -JVfPx/po - dU{x)/dx + Ex{t) + yj2^kT^x{t), 
Py = -IVFPy/PO + Ey{t) + ^2^kTiy{t), 
V = {x,y) = de/dp = vfp/po, (1) 

where 7 denotes a phenomenological constant damp- 
ing coefficient The random forces £,x{t) and (,y{t) 
are two white Gaussian noises with {Ci{t)) = and 
(^i(i)^j(O)) ~ Sij5{t), for I, J = x,y, which ensures proper 
thermalization at temperature T. We recall here that 
the modulus of V is constant and cc^ual to 0- These 
equations describe a relativistic Brownian dynamics, thus 
suggesting new settings for the experimental and theo- 
retical investigation of relativistic thermodynamics and 
relativistic Brownian motion. Below we report results 
from extensive numerical simulations we performed by 
integrating Eqs. ([T]). 

As anticipated above, the validity of Eqs. ([ij for 
electrons in graphene is restricted to the quasi-classical 
limit, where quantum mechanical effects can be safely 
neglected. Thus, the testing ground for relativistic phe- 
nomena on a graphene chip is subject to the following 
physical restrictions (isj : 

(i) the de Broglie wavelength across the barrier, A^; = 
h/px, must be much smaller than the period of substrate 
potential, A^; <^ L, so as to neglect either miniband or 
discrete energy levels; 

(ii) the de Broglie wavelength parallel to the barrier, 
A J, = h/py, must be much smaller than L, so as to sup- 
press the probability of chiral tunneling (Klein paradox 

The condition Xy <^ L allows neglecting, in par- 
ticular, the electron-hole tunneling; 

(iii) temperature (or noise strength) must be high enough 
to further ensure that chiral tunneling at the top of the 
U{x) barriers is negligible with respect to the competing 
noise-activated particle hopping, i.e., 

mm{hvF/L; hulU^ /vfA^L} < kT < U. 

Condition (iii) on T is introduced [l^l because noiseless 
carriers moving along the x-axis with Xy <^ L would still 
be subject to Klein's mechanism (barrier transparency) 
Note, however, that in 2D geometries with smooth 
potentials U{x), the fraction of electrons undergoing 
Klein's mechanism is negligibly small, not only at high 
temperatures [as guaranteed by condition (iii)], but also 
in the presence of sufficiently strong transverse drives. 
Indeed, drives oriented along the y axis tend to suppress 
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FIG. 2: (Color online) Relativistic ratcheting in an asym- 
metric, periodic potential U{x) = sin (27ra;) + 0.25 sin (47ra;), 
driven hy E = (0, Acos {ujyt)): (a) Net cmrent Jx/vf versus 
T for different LUy (reported next to each curve); (b) Rescaled 
diffusion coefficient, Dx/Do, versus T for the same simula- 
tion parameters as in (a). Here, Do = fcT/7 and AU ~ 2.20 
(height of the substrate barriers). Other parameters are 
A = 20 and vf = j = k = 1. 

chiral tunneling even for particles with py = [l5| , thus 
corroborating our quasiclassical model. 

Stochastic effective mass - The moment components 
Px and Py in Eqs. ([1]) are coupled via the dispersion 
relation, e = vppo, even in the absence of a substrate 
potential, U{x) = 0. This gives rise to a peculiar phe- 
nomenon, which characterizes the relativistic Brownian 
motion with respect to its non-rclativistic counterpart 
[Fig. Iljb)]. By numerically solving Eqs. ([T]) for electrons 
driven, say, along the y axis hy E = (0, Acos (ojyt)), we 
determined the strength, S, of the delta-like peak of the y 
power spectral density at the driving frequency, LOy. Con- 
trary to the non-relativistic limit, where it is T indepen- 
dent, at high temperatures S decreases according to the 
approximate power law T~^. This result can be quantita- 
tively explained by noticing that according to the high- 
T Jiittner distribution, the particle energy distribution 
is yo(e) ~ eexp(— e/fcT); hence, (e) = vf{po) = 2kT. 
Upon replacing po by 2kT/vF, Eq. ^ is reduced to a 
set of decoupled Langevin equations describing a non- 
relativistic particle of effective mass mofi = 2kT/vp. Ac- 
cordingly, 5/5*0 = + i^y/'^o)'^] with So = (A/^LUy)'^ 
and Wo = l/mcs = v\^l2kT [dashed line in Fig. [Hb)]. 
Stochastic mass renormalization provides a simple vali- 
dation check for our phenomenological model in Eqs. ([Ij. 

Relativistic ratcheting - We consider next charge car- 
rier transport due to the rectification of non-equilibrium 
perturbations on a substrate with asymmetric potential 
U{x) (ratchet or Brownian motor effect [iBl)- The work- 
ing principle of a Brownian motor is that, under certain 
conditions, asymmetric devices are capable of rectifying 
random (i.e. noisy) and/or deterministic (periodic) mod- 
ulations. In the non-relativistic regime, ratcheting occurs 
only in the x direction, as a mere effect of the x com- 



ponent of the periodic drive, Ex{t). The rectification 
current weakens if E{t) is rotated at an angle with the 
X-axis, until it drops to zero for ac drives parallel to the 
substrate valleys. Indeed, the component Ey(t) of the ac- 
field does keep the system out of equilibrium, but cannot 
be rectified, since the substrate potential is uniform in 
the y direction. Stated otherwise, the x and y dynamics 
are decoupled. In the relativistic Langevin equations ([T|), 
instead, the orthogonal ac-drive components, Ex{t) and 
Ey it) , are nonlinearly coupled through the dispersion re- 
lation, so that both can be rectified by the asymmetric 
potential U{x). Most remarkably, apart from a special 
parameter range dominated by relativistic chaotic effects 
(see below), transverse rectification induced by Eyit) in 
the X direction only occurs at finite temperatures, thus 
implying a bona fide noise-sustained Brownian ratchet 
[la ]. Noise is required to force particle fiuctuations in 
the x-direction, around the asymmetric minima of U{x); 
as Py is driven by Ey{t) toward values of the order of 
Vt, or smaller, correspondingly Vx jumps to a maximum, 
Vx ~ vfPx/\Px\ = ivp, so that the particle is kicked in 
the X direction, either to the right or to the left; this is 
how spatial asymmetry comes into play [l7j . Naturally, 
such a rectification mechanism becomes ineffective when 
the particle sits at a potential minimum with px = 0, 
which only occurs in the absence of noise, i.e., for T = 0. 

To demonstrate the relativistic transverse ratchet ef- 
fect, it suffices to arrange the electrodes on the graphene 
chip so as to generate an asymmetric periodic potential, 
U (x) , directed along x and ac drive the charge carriers 
along y with E{t) = {0, Acosujyt). By numerically inte- 
grating the corresponding Eqs. ([T]), we observed that a 
relativistic particle tends, indeed, to drift in the x direc- 
tion with net current Jx = (x) and diffusion coefficient 
Dx = \imt^ao[{xitf) - {x{t))'^]/2t. The dependence of 
Jx and Dx on the temperature and the driving frequency 
are displayed in Fig. [5J In panel (a) it is apparent that 
at constant ujy the absolute value of the current, | |, 
hits a maximum for an optimal T and tends to zero in 
the limits T — >■ and T — >■ 00, thus underscoring the key 
role played by noise. This result supports our conclusion 
that we deal with a new category of noise-sustained (no 
current at zero noise) and relativistic (no ratchet phe- 
nomenon in the non-relativistic limit) Brownian motors. 
The direction of the net rectification current depends on 
both the driving frequency and the temperature; its sign 
can be (multiply) reversed by changing either T, panel 
(a), or ujy (not shown). Most remarkably, in panel (b) 
one can locate a finite temperature interval, where, for a 
given cjy, the scaled diffusion coefficient Dx/Dq decreases 
with increasing T; indeed, Dx approaches the expected 
asymptotic value Dq = kT/j only after going through a 
minimum. This T-interval approximately coincides with 
the T-interval in panel (a), where the absolute value of 
the current and, therefore, the mobility of the charge car- 



4 




FIG. 3: (Color online) Relativistic stochastic resonance, (a) The confining potential, U{x) — x^/A — x"^ /2, bistable in the 
i-direction and uniform in the transverse direction y\ (b) Corresponding phase diagram for E = {Q, Acos{u)yt)). Insets: 
trajectory samples for points in the plane (^,0;^) (denoted by arrows) belonging to different dynamical phases, (c) Residence 
time distributions, p(r), for different T (as indicated) and constant drive, A = 1^ and ujy — 0.5. Relevant simulation parameters: 
integration step, f.2 x fO^^; run time length, 5 x fO^; distribution time bin, 2.5 x 10~^. Inset: the distribution peak strength, 
S = J^^ p{T)dT, plotted versus T (SR signature). Other simulation parameters are: vf ~ "f = k = 1. 



riers is the largest. Analogously, we checked that the ab- 
solute maximum of as a function of ujy corresponds to 
a local minimum of Jx versus ujy. This observation sug- 
gests that for appropriate system parameters, elementary 
charge excitations in graphene can be regarded as Brow- 
nian particles with large mobility and low diffusivity. 

Chaos and stochastic resonance - To investigate the 
role of anomalous signal amplification we consider rel- 
ativistic Brownian motion driven along y by the time- 
periodic electric field E{t) = (0, A cos{ujyt)), and confined 
along X by the double- well potential U{x) = /A — x'^ /2, 
sketched in Fig.[3l^a). In the deterministic, or noiseless 
regime, T = 0, the motion of a classical particle in the 
x-dircction would be frozen, i.e. the particle would sit 
in one of the two U {x) minima, regardless of the drive 
applied in the y direction. In clear contrast, the re- 
sponse of the relativistic particle in the z-dircction is 
strongly affected by the nonlinear coupling of the a;- 
and y-dynamics, controlled, respectively, by U{x) and 
Ey{t). In fact, by appropriately tuning the drive pa- 
rameters A and ujy, the relativistic particle can execute 
either small amplitude intra-well oscillations within one 
confining well, or large-amplitude inter-well oscillations 
between the two confining wells, or even exhibit chaotic- 
like switching between the two wells. The phase diagram 
of these three different dynamical regimes is depicted in 
Fig.El^b) with samples of the corresponding x-trajectories 
reported in the insets. 

Our numerical simulations showed that in the chaotic 
regime a relativistic particle performs small amplitude os- 
cillations around one U (x) minimum for an unpredictable 
time interval (residence time) before switching over into 
the other minimum, where it resumes its small ampli- 
tude oscillations until the next switching event. These 
switchings are due to the intrinsic instability of the par- 
ticle dynamics and not to thermal fluctuations, as in our 



simulations the temperature was initially set to zero. 

On raising the temperature, the inter-well dynamics 
induced by the transverse field, Ey{t), becomes increas- 
ingly noise dominated. In Fig.|3ljc) we display the nor- 
malized distribution density of the residence times in ei- 
ther U{x) well, T, at different T. The peak split structure, 
characteristic of the noiseless chaotic dynamics, merges 
into a single broad peak centered at around r = 2-k /ujy. 
The corresponding peak strength, S*, defined here as the 
area encircled by the peak, attains a maximum for kT 
smaller than, but close to the barrier height of U{x)^ 
AJ7 = 0.25. Such a behavior can be regarded as an 
instance of stochastic resonance, a phenomenon well es- 
tablished for bistable non-rclativistic Brownian motion 
[lit , with two additional peculiarities: Here the periodic 
drive acts perpendicularly to the hopping direction and 
the resonance is related with chaotic rather than periodic 
switches between potential wells. 

We conclude by remarking that, irrespective of the 
phenomenological details of Eqs. ([T]), the spatial coor- 
dinates of a relativistic particle in 2D are nonlinearly 
coupled via a non-separable energy-momentum disper- 
sion relation. In this regard our approach may apply not 
only to graphene but also to other materials, like certain 
types of semiconductors. Moreover, the results discussed 
in this work, besides suggesting an inexpensive technique 
to assess the current stochastic formalisms for relativis- 
tic particles, might well find application in ultrafast elec- 
tronics. Indeed, pursuing the technological implemen- 
tation of the effects reported above (and yet others as 
summarized in Rcf. jl3l] ) is expected to become viable, 
should one eventually succeed in replacing GaAs- with 
graphcne-based electronics [l^l- 
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o : 

^ ^ , I. QUASI-CLASSICAL DYNAMICS OF CHARGE CARRIERS IN GRAPHENE 

X) . 

D ■ . 

I ■ We discuss here our arguments in support of the semiclassical description, Eqs. (1), main text, of the carrier 
' dynamics in a reahstic graphene sheet coupled to a heat-bath at given temperature T. 

, In the low-energy limit and in the absence of randomness (fluctuations, T = 0, and disorder, alike), the charge 
carriers in graphene, electrons and holes, obey the Dirac equations, 

^ , ■' VF VF Ot 

{Vx + iVyj'i/A H Wb = ?r— , 

I , vf vf at 



where the two spinor components of the carrier wave function, refer to the two triangular graphene sub- 

^ \ lattices and the momentum operators are defined as {VxjVy) — {—ihd/dx, —ihd/dy). 

Applying the operators (Vx — iVy) and (Vx + iVy), respectively, to the first and the second equation leads to the 
coupled second-order partial differential equations, 



o 



F . 

f^(^ ^\ 2—U—- — 

\ 9x2 Qy2 J y2 Q-f-2 y2 y2 



Vp dt Vf \dx dy 
ih dU ih ^ f dU dU 
Vp dt Vf \dx dy 



Equivalently, the effective actions [SatSb), defined by rewriting a,'^ b) as (exp(iS'^/?i), exp(zS'B/?i)), satisfy the 
coupled equations. 



> 

■ \ dx J \ dy J Vp\dt J dx^ dy'^ Vp dt^ Vp dt vf \ dx dy 
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ih d^SA 
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^ vl dt^ 
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d^SB 


^^d'SB 


ih d^SB 


ih dU 


dx^ 


dy"^ 


^ vl dt^ 


'^~dt 



, \ dx J \ dy J Vp \ dt J dx^ dy^ Vp dt^ Vp dt vf \ dx dy 



In the limit ?i — > 0, Sa and 5b tend to the unique classical action, S, solution of the classical Hamilton- Jacobi 
equation, 

dsV fdsV 1 fdS 



On making use of the well-known relation, % — dS/dt, one eventually derives from Eq. the classical Hamiltonian 
for a relativistic massless particle. 



n = ±VF^Pl+Pl+u, (4) 

with U{x,y,t) = U{x) — Ey{t)y — Ex{t)x and ± referring respectively to electrons and holes, whose densities are 
related by electroneutrality. 

Note that the lowest order corrections to the classical limit, Eq. ([3]), are 0{h), at variance with the standard WKB 
formalism for the Schrodinger equation. Such corrections, for \Sa — Sb \ ~ h, are due to the electron- hole transitions, 
so that a systematic quasiclassical treatment would require a two particle species description. 
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Thus far, the interactions of a single charge carrier with other carriers, phonons and impurities have been neglected. 
If the relevant collision times are short compared with the operation time constants of the graphene chip (periods of 
the ac-drives, carriers drift times across the substrate potential cells, etc.), then such interactions can be accounted for 
by adding to the potential function, U, an appropriate fluctuating term /7^, namely, U{x,y,t) — U{x) — Ey{t)y — 
Ex{t)x + U^, which corresponds to an isotropic white noise force in the plane {x,y). 

The same interactions are also responsible for the damping term in the Langevin equations (Eq. (1), main text). 
The appropriate form of the damping and the fluctuating terms are determined by the fluctuation-dissipation theorem. 
In Eqs. (1), main text, we first introduced a damping term of the viscous type, —^V [2], and then determined the 
random forces accordingly Q. Such a choice is in line with the single r-approximation of the Boltzmann collision 
integral, which, for systems not too far from equilibrium, naturally yields a viscous term with 7 = 1/r, see, e.g., in 
Ref. [4]. 

The damping coefhcient 7 contains contributions from all possible collision mechanisms a charge carrier may undergo 
when moving through a graphene chip. To determine the relative importance of such mechanisms would require a 
detailed quantum kinetic theory (see, e.g. Ref. [1]), which rests beyond the purposes of our semi-classical approach, 
and, more importantly, direct access to experimental parameters, such as spatial disorder (impurities), internal lattice 
stress in the presence of a substrate (phonon distribution), local temperature and magnetic field fluctuations, which 
are hardly controllable. Recent experimental evidence 0, 0] seems to indicate that the main source of carrier viscosity 
in graphene is attributable to electron-phonon interactions. Anyway, at the present stage, 7 ought to be regarded as 
a phenomenological parameter of the semiclassical Eq. (1), main text, to be experimentally determined. 

The semiclassical approach outlined above provides a meaningful phenomenological description of charge carrier 
transport in a graphene chip as long as its characteristic length scales (e.g., the chip size and the electrode spacings) 
are much larger than the de Broglie wavelengths, A^^ = h/p^ and Ay = h/py, with {px,Py) — {dS/dx,dS/dy), see 
conditions (i) and (ii) in the main text. Condition (i) states that for A^; much shorter than the characteristic length 
scale of the substrate potential along the x-axis, L ~ \U/{dU/dx)\ for smooth functions U{x), the discreteness of the 
carrier energy spectrum becomes negligible, as to be expected in view of the Bohr-Sommerfeld quantization rule. 

Condition (ii), h/py <C L, follows from requiring that the probability of chiral tunneling is distinctly smaller than 
unity, i.e., that the particle must be reflected many times by the potential barrier before it may tunnel across. Chiral 
tunneling is exponentially suppressed for non-zero transverse momenta, Py ^ 0, and smooth time-independent poten- 
tials, U{x). In this case, the classical action has the form S — Pyy — Et + So{x), where E is the conserved energy and 

So{x) can be readily derived from the the Hamilton- Jacobi Eq. ([3]), that is, Sq — ± J ^ {E ~ U)^/vp — Pydx. Chiral 

tunneling occurs through the classically forbidden region, \E — U\ < vplPyl; therefore, the smoother the potential, 
the larger the size of the forbidden region, 2vF\Py\/\dU /dx\. As a consequence, chiral tunneling is characterized by an 
exponentially small probability, w = ex'p{—(oVFPy/JT'\dU/dx\) @, with ^0 denoting a constant of the order of unity, 
which may be appreciable only for a small fraction of electrons with py ^ y^h,\dU/dx\/vF- By introducing the estimate 
Py ^ U/vp for the transverse momentum, we derived condition (ii). It is important to stress that a time dependent 
electric field acting along the barrier, Ey{t), tends to exponentially suppress also the chiral tunneling of electrons 
moving perpendicularly to the barrier, py = 0, that is, w — exp[— 7re^Wi?A^/(4?i|dJ7/da;|a;^)] [9^]. This remark further 
corroborates our quasiclassical approach in the presence of strong enough drives. 

As an additional consistency requirement for our semiclassical description of graphene charge carrier transport over 
the substrate potential, we must further ensure that chiral barrier tunneling is negligible with respect to the thermal 
barrier hopping, that is, w <^ exp{—U/kT). This yields the lower bound for temperature in condition (iii). 

In conclusion, our quasiclassical Langevin equations, Eqs. (1), main text, represent a valid starting point to 
model carrier transport in graphene at finite temperature (see e.g., jTo|). with one proviso: Such a phenomenological 
description only applies as long as (i) chiral tunneling (which implies the conservation of the associated chiral quantum 
number), (ii) interband transitions, and (iii) intervalley scattering can safely be ignored. As mentioned above, such 
effects could be accounted for, at least in principle, by developing a more refined phenomenological (semiclassical) 
model, which involves two (or more) species [III] of interacting particles, electrons and holes, and thus is capable to 
describe an even richer relativistic stochastic dynamics. 

II. RELATIVISTIC VS. NON- RELATIVISTIC BROWNIAN MOTION 

The new effects observed so far are summarized in TablelH where for reader's convenience we compare case by 
case the different behaviors of relativistic (massless) and non-relativistic (massive, or overdamped) Brownian particles 
confined to a two dimensional geometry. One has to bear in mind that the orthogonal coordinates of the relativistic 
particle described by the Eqs. (1), main text, are dynamically coupled even in the absence of a nonlinear substrate 
or external gradients. Such a coupling vanishes in the non-relativistic limit. As a consequence, first, the response 
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of graphene charge carriers to a time-periodic drive is temperature dependent, see Fig. lb, whereas the response of a 
classical free particle only depends on the drive. Second, on applying a periodic symmetric spatial modulation in a 
given direction, a relativistic massless particle can be rectified in that direction, when driven by two center-symmetric 
unbiased time periodic fields (harmonic mixing), acting in orthogonal directions. Note that our data on harmonic 
mixing are not reported in the manuscript. Rectification occurs only for a certain combination of the relevant driving 
frequencies In the case of a classical particle, harmonic mixing would require that both harmonic signals act 

in the direction of the spatial modulation. Third, in the presence of a spatially asymmetric ratchet-like potential, 
U{x) (Fig. 2), the relativistic particle can be rectified along x by periodically forcing it along y, contrary to classical 
ratchets, where rectification of transverse signals was never observed. Finally, if the motion of a noiseless relativistic 
particle is confined to a double-well potential (Fig. 3) along x, in the presence of a transverse harmonic drive, Ey{t), 
its coordinate x{t) can be either constant, or periodic, or chaotic, depending on the Ey{t) parameters. Adding noise 
causes an unusual manifestation of stochastic resonance with no counterpart in the non-relativistic regime: at finite 
temperature, particle switchings between U{x) minima can be optimally synchronized even by means of a transverse 
periodic modulation. 

This work was partly supported by the HPC-Europa2 Transnational Access Programme, application N 278; Euro- 
pean Commission (NanoPower project); the Alexander von Humboldt Foundation, EPSRC (EP/D072581/1) and the 
German excellence cluster" Nanosystems Initiative Munich" (NIM) via its seed funding program. 



TABLE I: Driven stochastic dynamics of relativistic and non-relativistic particles: A comparison 





Relativistic massless particles 


Classical massive or 
overdamped particles 


Experimental data sources 


Cosmic rays, graphene chip 


nano-particles in liquids 
and nano-channels, 
vortices in superconductors, 
bio-molecules in bio-systems etc. 


Possible applications 


THz graphene-based electronics 


nanotechnology, 
noise-driven nano-robots, 
drug delivery etc. 


Current state of the art 


conflicting formulations 
limited applicability 


well established working tool 


Coupling between orthogonal coordinates 


Yes 


No 


Rectification in a spatially symmetric 
periodic potential U{x), induced 
by mixing of two harmonic signals 


Occurs if the signals are 
applied in the same as well as 
in the orthogonal directions 


Occurs only if the two signals are 
applied in the same direction 


Rectification in a spatially asymmetric 

periodic potential U{x) induced by 

a single center-symmetric unbiased signal 


Occurs if the signal is applied along 
the X, as well as along the y direction 


Occurs only if the signal 

is applied along the x direction 


Motion in a double-well potential 
with orthogonal harmonic drive; 
chaos and resonance 


Zero temperature lateral response 
can be constant, periodic or chaotic 
in time. 

Finite temperature response in the 
chaotic regime: non-exponential 
double-peak distribution of the 
residence times; stochastic resonance 


No zero temperature lateral response 

Finite temperature response: random 
switchings with exponential distribution 
of the residence times 
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